Oscillon resonances and creation of kinks in particle collisions 
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We present a numerical study of the process of production of kink-antikink pairs in the collision 
of particle-like states in the one-dimensional 4 model. It is shown that there are 3 steps in the 
process, the first step is to excite the oscillon intermediate state in the particle collision, the second 
step is a resonance excitation of the oscillon by the incoming perturbations, and finally, the soliton- 
antisoliton pair can be created from the resonantly excited oscillon. It is shown that the process 
depends fractally on the amplitude of the perturbations and the wave number of the perturbation. 
We also present the effective collective coordinate model for this process. 
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Introduction. Non-linear field theories in the weak 
coupling regime usually contain two different mass scales 
associated with the perturbative particle-like states and 
with the soliton sector of the model, respectively. Conjec- 
ture about the role nonperturbative effects, related with 
the production of the soliton-like states in particle colli- 
sions, may play in high energy physics has attracted 
a lot of attention recently. Over the last two decades the 
problem of the transition between perturbative and non- 
perturbative sectors of the theory has been considered in 
several contexts. 

Simplest example of the topological solitons in one di- 
mension is the kink solution of the 4> A model. This model 
has a number of application in condensed matter physics 
0, field theory 0,|| and cosmology [|[ . Dynamical prop- 
erties of kinks, the processes of their scattering, radiation 
and annihilation have already been discussed in a number 
of papers, see e.g. [§-[l6j]. In integrable theories, like the 
sine-Gordon model, there is no energy loss to radiation 
and kinks do not annihilate antikinks. However in the 
non-integrable 4 model, the radiation effects in the pro- 
cess of kink-antikink {KK) collision become very impor- 
tant and depending on the impact velocity, the collision 
may produce various results, e.g., an oscillating bound 
state can be formed, also the soliton and antisoliton may 
bounce and reflect from each other. 

Although the process of annihilation of the KK state 
of the (p 4 model has been investigated in detail [Til . H~3l . 
[IH , there is not much information about the inverse pro- 
cess, the creation of the KK pairs by the collision of 
two identical bunches of particles. In the recent work § 
production KK pair was considered in assumption that 
two colliding wave trains are composed of the bunches of 
identical breathers, i.e., tightly coupled KK states. Ev- 
idently, the kink-antikink production may proceed even 
in the case when there is no kink-like states in the ini- 
tial configuration at all. Here we aim to elucidate the 
mechanism for this process. 

In is known that the collision of a kink and an antikink 
is chaotic, i.e., for some values of the impact velocity 
the solitons bounce back while for some different impact 
velocity, smaller or larger, they annihilate [J, [16|. This 
behavior is related with a resonance effect between the 



oscillations of the KK pair and excitation of the discrete 
vibrational mode of the kink. 

So we might expect the opposite process of the pro- 
duction of the KK pairs in the collision of particles also 
will have similar fractal character due to resonance effect 
between the oscillon created in the particle collision and 
the oscillation of the correlated KK pair. 

Here we investigate the oscillon resonance numerically. 
We find that this resonance excitation plays a crucial role 
during the process of creation of the KK pairs in the 
collision of particles. We observe furthermore the fractal 
structure of this process. 

The model. We consider the standard one-dimensional 
(j) 4 theory, with two vacua <f> = ±1, defined by the rescaled 
Lagrangian density 
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The perturbative sector of the model consists of small 
linear perturbations around one of the vacua with the 
mass m = 2. The static kink solution for this model 
interpolates between the vacua 4>q — — 1 and 4>q = 1 as x 
increases from — oo to oo: <j}x{x,t) — tanhx. 

There is a lot of similarity between the the non- 
integrable model ((I) and its integrable sine-Gordon coun- 
terpart. However, the states of the perturbative sector 
are different in these theories. Evidently, in both mod- 
els there are zero translational modes in the spectrum 
of the linear perturbation about the kinks but a single 
<j) kink has in addition a normalizable discrete vibra- 
tional mode which oscillates harmonically with frequency 
u)\ = The continuum modes on the kink background 
have higher frequencies uj > 2. Evidently, if the ampli- 
tude of the oscillation is large enough, such a periodi- 
cally expanding and contracting kink can be treated as 
kink-antikink-kink bound state and this excitation can 
be considered as an intermediate step in the process of 
creation of the KK pair on the kink background [HI, [ll| . 

Another situation is related to the possibility of pro- 
duction of the KK pairs on the trivial background. In- 
deed, the linear excitation spectrum around the trivial 
vacuum contains the radiation modes and within the </> 4 
model, the collision of these particle-like states may pro- 
duce KK pairs Q. 
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Note that non-linear field theories usually contain sev- 
eral types of topological and non-topological excitations. 
Indeed, besides the solitonic configurations there is an- 
other spatially localized non-perturbative oscillon solu- 
tion which, although unstable, are extremely long-lived 
[l8l - l20l |. The oscillon states naturally appear in various 
models (H^. 

In the 4> 4 model the oscillon solutions are almost pe- 
riodic. One can find the oscillon numerically by solving 
the field equation in the Fourier series in time: 

cf>=l + rj (x) + rn(x) cos(fti) + 772(2) cos(2fit) + . . . (2) 

If 57 < m — 2, the oscillations are below the threshold 
and cannot propagate as modes of the continuum, so the 
oscillon remains relatively stable and the r\\ term domi- 
nates. 

It was pointed out recently that an oscillation mode of 
the 4 model may decay into a KK pair Q. 

Numerical results. The initial data used in our sim- 
ulations represent two widely separated identical wave 
trains propagating from both sides on the trivial back- 
ground towards a collision point: 

4>(x, t) = l+C[F(x+vt) sin(u)t+kx)+F(x-vt) sin(ujt-kx)] 

where k is the wavenumber of the incoming wave, w = 
\/k 2 + 4 is the frequency and v = k/ur is the velocity 
of propagation of the wave train. We consider the en- 
velop of the train F(x) — [tanh(a; — a\) — tanh(x — 02)], 
also the gaussian envelope F(x) — e - a ( x - a ^) wa s used 
to prove that our results are independent of the particu- 
lar choice of the initial state. The parameters ai,<Z2 and 
03 define the length of the train and the initial separa- 
tion between the trains. Typically, we used the values 
a\ = 10,a2 = 30,a3 = 20. The amplitude C and the 
wavenumber k are the impact parameters, which can be 
changed freely. In our numerical analysis we found that 
after small amplitude collisions, the two wave trains sep- 
arate and move in opposite directions and the radiation 
is created due to the interaction between these trains. In 
the center of collision an oscillating lump remains. For 
small amplitudes, the frequency of the oscillation is just 
a bit above the mass threshold. This indicate that the 
lump could be identified with low wave-number linear 
excitation of the trivial vacuum. For large amplitude 
collisions, the remaining lump oscillates with frequency 
within the mass gap, so such a state can be identified as 
an oscillon. 

Furthermore, for a certain range of values of the impact 
parameters, C and k, we observed creation of the KK 
pairs. During this process also an oscillon is create in 
the collision center (Fig. [TJ. The most important feature 
of this process is that in the space of parameters, the re- 
gions of creation of the solitons and the regions where this 
process in not taking place, are separated by a fractal- 
like boundary (Fig. [2]). For finite wave trains we do not 
expect this boundary would be a real fractal but some 
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Figure 1. Production of the kinks in the collision of two iden- 
tical wave trains. The initial and final field configurations are 
plotted at t — 0, t — 45 and t = 100 respectively. 
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Figure 2. Fractal structure in the C, k plane. Shading (or 
colour) represent the measured minimum of average of the 
field (A) = f^_ 10 dx <fi(x,t). The dark regions (blue in 
colour), where (A) < — 1 indicate creation of the KK pairs. 

properties of scaling are observed (Fig. We have mea- 
sured the fractal (box) dimension to be d — f .770±0.0f 1, 
which is much more than 1. In the case of the gaussian 
envelope we found d = 1.865 ± 0.007. The interesting 
peculiarity of the latter case is that the KK pairs can be 
created even if k = (standing wave perturbation). 

For certain values of impact parameters, an oscillon re- 
maining in the collision center decays into the second KK 
pair. Sometimes the second pair moves even faster than 
the first pair, and it may annihilate with the first one 
creating two moving oscillons. We know that the process 
of collision KK pair also leads to fractal structure in the 
velocity space [lj],[l6j]. ^ n our process, instead of creation 
on KK pair, two oscillons could also be ejected from the 
collision center and after a while they could decay into 
two pairs of KK. We observed some evidences that these 
processes also yield the fractal dependency of impact pa- 
rameters. This study will be reported elsewhere. 

To find a numerical solution of the PDE describing 
the evolution of the system, we used the pseudo-spectral 
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The Lagrangian of the free oscillon has the form 



Figure 3. Calculation of the box fractal dimension. Number 
of the boxes N covering the boundary is plotted versus the 
size of the covering boxes I. The slope in the above log — log 
plot gives the dimension d = 1.770 ±0.011. 



method. For the time stepping function we used sym- 
plectic (or geometric) integrator of 4th order to ensure 
that the energy is conserved. 

Effective collective coordinate model. In order to 
capture the most important steps in the process of the 
creation of KK pair, in the collision of two identical 
bunches of particles, we use the collective coordinate 
method which allows us to identify the physical degrees 
of freedom of the system under consideration. This ap- 
proach has been applied to describe the dynamics of the 
kink-antikink system [17] . 

First, we describe the process of creation of the oscillon 
in the collision of the incoming wave trains. We assume 
an initial field configuration on the trivial background 



(f>(x, t) = 1 + A(t)/ cosh(x/a;o) + £0, t), 



(4) 



which corresponds to the profile of the oscillon solution 
[l3l [l9| with some additional perturbation £(x,t). The 
gaussian approximation to the oscillon configuration [2(| 
was also used to check the results. Here the variable 
A(t) is introduced as the collective coordinate of the os- 
cillon and the parameter xq represents the oscillon width. 
From the expansion (J2J we know that when £ = the 
oscillon should, in the first approximation, oscillate as 
A(t) = Aocos(ftt), where Aq is the amplitude of the 
oscillations, < 2 and the value of the parameter Xq 
depends on the amplitude Aq. In the presence of the 
external field £, the amplitude of the oscillon changes. 
However, for the sake of simplicity we set xq = 1.5 as 
it is the width of the oscillon oscillating with amplitude 
Aq = 0.4. Substituting (TJ| into ([T]) and after integration 
over all space x gives the effective Lagrangian which can 
be split into three parts: Lagrangian of the free oscil- 
lon, Lagrangian of the perturbation £ and Lagrangian of 
interaction between the oscillon and the perturbation, 
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This is the Lagrangian of an anharmonic oscillator with 
frequency = ,/4+ > 2. Since the frequency of 

V Jx 

the oscillon must be smaller than m — 2, the amplitude 
of the oscillations must be large enough to decrease the 
oscillation frequency below the mass threshold [2(J Hl[ , 
so the nonlinearities are crucial for the existence of the 
oscillon. 

We assume that the field £, is a solution to the equation of 
motion of the Lagrangian Lj. The perturbation £, should 
represent two wave trains coming from ±oo. For the sake 
of simplicity we take the perturbation of the form (J3J . 
In the last part of the Lagrangian we take only the linear 
term in £: 

L mt = a(t)A(t) + p{t)A 2 + l(t)A + 0(A 3 ) + 0(£ 2 ). (7) 



where 



a(t) = dx 
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£ x smh(x/x ) 
cosh(x/x ) ' x cosh 2 (x/a; ) 
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and 



6£ dx 



cosh (x/xq) ' 



6 dx 



cosh(x/a;o) ' 



(8) 



(9) 



L(A, A) = L A + L 6 + L 

i 



(5) 



although we have also investigated the effect of the higher 
order terms. Because £ is an oscillating function over a 
compact support, the above Lagrangian along with La 
describes dynamics which is similar to the dynamics gov- 
erned by Mathieu's equation, (especially the term pro- 
portional to /3(f)) but with additional nonlinear terms 
and a source term (proportional to a(t) and j(t)). We 
have studied the dynamics of the system numerically. 
Unfortunately we could not find the analytic form of 
the integrals a(t), f3(t) and 7(i), so we used numerical 
methods which are on the same level of complexity as 
the explicit solution of the underlying PDE. However, 
within the collective coordinate approach we can sepa- 
rate the most important degrees of freedom and show 
that the corresponding nonlinear interaction is responsi- 
ble for generation of the fractal structure. 

The initial condition is that A(0)=0. As the wave train 
approaches the point of the collision, the oscillon mode is 
excited. If the amplitude of the perturbation is relatively 
small, then the oscillon, created in the collision, oscillates 
with a constant amplitude around A = 0. However, if 
the amplitude is large enough, or the incoming pertur- 
bations are close to one of the (Mathieu) resonances, the 
amplitude of the oscillon rapidly increases and it starts 
to oscillate around A = — 1 (or, in other words, around 
4> = 0) with amplitude of order 1, as illustrated in Fig 
[U This clearly breaks our effective approach, but it also 
means that the system has changed the ground state. 
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Figure 4. Resonance of the amplitude of the oscillon (dashed 
line) within the effective model. 

Such a resonant oscillation with a large amplitude, on 
the other hand, shift the center of the oscillation. This 
transition can be related with creation of the KK pairs 
although the corresponding collective coordinates are not 
presented in our simple model ([5]). Again, when we ex- 
amined this effective model, we found a fractal structure 
on the plane A, k. This fractal structure was less com- 
plicated and more localised than in the case of full PDE. 
That means that although our effective model works and 
captures qualitatively the most important features of the 
full system it also fails to reproduce some of the details, 
which is not a surprise for such complicated dynamical 
process. We have also introduced an approximation for 
a(t), (3{t) and ^{t) and again, we could reproduce both 
the resonance excitation of the oscillon and the genera- 



tion of fractal structure. 

This result shows, that even after performing so many 
simplifications we could reproduce (at least qualitatively) 
the most important features of the evolution of the sys- 
tem. This result confirms our conjecture about the mech- 
anism of the creation of the KK pair in a 3-stage process 
(i.e., excitation of the oscillon, resonance and oscillon de- 
cay into the KK pair). Secondly, we conclude that the 
interaction between the incoming wave trains and the os- 
cillon is the underlying reason for the generation of the 
fractal structure. Thirdly, given the generality of our 
approach, we expect that the effective nonlinear interac- 
tions of the same type can be found in many different 
models, so the fractal structure should not be limited 
only to the case of the <p A model. 

Concluding, we have found that the KK pairs can be 
created from the collision of particles in the 3-stage pro- 
cess. In the first stage, the collision process produces an 
oscillon excitation. Next, the oscillon interacts with the 
incoming trains, and due to the parametric resonance it 
may decay into the KK pair. This process yields the 
fractal structure in the space of parameters of the model. 
We reproduced the most important features of the evo- 
lution of the system by constructing an effective model. 
Our numerical investigations indicate fairly clearly that 
the observed phenomena seems to be independent of the 
approximations we used. We expect that other topolog- 
ical defects, like vortices, hopfions or monopoles, can be 
created in the process of the same type, i.e., via resonance 
excitation of an oscillon created by the particle collision. 
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